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Abstract 

A new class of lattice Dirac operators which satisfy the index theorem have 
been recently proposed on the basis of the algebraic relation 75(75.0) + (75/5)75 = 
2a 2k+l (^D) 2k+2 . Here k stands for a non-negative integer and k = corresponds 
to the ordinary Ginsparg- Wilson relation. We analyze the chiral anomaly and index 
theorem for all these Dirac operators in an explicit elementary manner. We show 
■ that the coefficient of anomaly is independent of a small variation in the parameters 

r and mo, which characterize these Dirac operators, and the correct chiral anomaly 
£^ ' is obtained in the (naive) continuum limit a —* 0. 

o ' 
in . 

1 Introduction 

o 

A new class of lattice Dirac operators D have been recently proposed on the basis of the 
algebraic relation [1] 

Oh" 75(75^) + (75^)75 = 2a 2k+ \ 1 ,D) 2k+2 (1.1) 

<D ! 

where k stands for a non-negative integer, and k = corresponds to the ordinary Ginsparg- 
Wilson relation [2] for which an explicit example of the operator free of species doubling 
has been given by Neuberger[3]. It has been shown in [1] that we can in fact construct 
the lattice Dirac operator, which is free of species doublers, for all values of k. Here 75 is 



a hermitian chiral Dirac matrix and 75-D is also hermitian. 
When one defines 

r 5 ee 75 - (a l5 D) 2k+1 (1.2) 

the relation (1.1) is written as 

r 5 (7 5 £) + (7 5 £)r 5 = 0. (1.3) 

The index relation [4] [5] on the lattice is generally written as 

T7T5 = n + — n_, (1.4) 

which is confirmed by[l] 

T7T5 EE J^^^n 
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A„=0 A n ^0 

= E ^n r 50n 
A„=0 

= E A[75 - (a 75J D) 2fc+1 ]0n 

A„=0 

= n + — n_ = index (1.5) 
where n± stand for the number of normalizable zero modes in 

75^0n = (1.6) 

for the hermitian operator 75.D with simultaneous eigenvalues 75 <p n = ±0„. We also used 
the relation following from (1.3) 

7 5 7jr 5 n = -A„r 5 n (1.7) 

if 

75 7J0 n = A n n , (1.8) 

which suggests that either T 5 (f> n for A n 7^ is orthogonal to <p n or else T 5 (f> n = 0. The 
positive definite inner product is defined by summing over all the lattice points 

0n0n = (0n> <t>n) = E ^ 4>l{x)<Pn{x) (1.9) 

x 

but the coordinate x is often omitted in writing <f> n . 

The Euclidean path integral for a fermion is defined by 

J VtpVijexplj i)Di)] (1.10) 

where 

/ ^D^ = ^(x)D(x,y)^(y) (1.11) 

and the summation runs over all the points on the lattice. The relation (1.3) is re-written 

as 

75 r 575 /j + /jr 5 = o (i.i2) 

and thus the Euclidean action is invariant under the global "chiral" transformation [5] 

4>(x) -> 4>'(x) = 4>(x) + i J24>(z)tl5F 5 (z, x) 75 

z 

Hy) - 1>'{y) = 1>(y) + * E er s(i/ ? «#H (i-is) 

w 

with an infinitesimal constant parameter e. Under this transformation, one obtains a 
Jacobian factor 

v^'vi)' = jvi>Vtp (i.i4) 
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with 

J = exp[— 2iTreT 5 ] = exp[— 2ie(n + — n_)] (1-15) 

where we used the index relation (1.5). This derivation may be regarded as a lattice 
counter part of the continuum path integral [6]. 

In Ref.fl] it was shown by using the method in [7], which is a lattice extension of 
the method in [6], that the index n + — n_ appearing in the Jacobian factor is related 
to the Pontryagin number for any operator in (1.1) if the operator 75.D satisfies suitable 
conditions. In this paper, we evaluate T7T5 in a more explicit and elementary manner 
on the basis of explicit formulas for 75 D in the continuum limitQ. We show that these 
operators 75-D for all k in fact reproduce the correct chiral anomaly and consequently 
correct Pontryagin number. 



2 A brief summary of the model and notation 

The operator Ts appearing in the index relation (1.5) has an explicit expressionfl] 

T 5 = 7 5 - F( 2 fe+i) (2- 1 ) 

with 

H (2k+1) = h 5 aD) 2k+1 = - 76 [1 + D w k+1) ^= 1 1. (2.2) 

y{L> w pv w 

The operator D w k+1 ^ is in turn expressed as a generalization of the ordinary Wilson Dirac 
operator as 

D % k+l) = i{gy k+1 + {B) 2k+1 - (2.3) 

See Appendix for further details of the general solution to (1.1). 

The ordinary Wilson Dirac operator Dw, which corresponds to Dy? •, is given by 

D w (x, y) = ij^C^x, y) + B(x, y) - -m Q 5 Xty , 

a 

C^(x,y) = —[Sx+fia.yUM - 5 X}y+(ia Ul(x)], 

B (x,y) = 7rJ2i 2S ^y ~ 5 V+^,xUl(x) - Syje+fiaU^y)], 

U^{y) = expliagA^y)], (2.4) 

where we added a constant mass term to Dw- Our matrix convention is that 7^ are 
anti-hermitian, (7^)^ = — 7 M , and thus (? = ^C^in^rn) is hermitian 

<?t =q. (2.5) 

Since the operators (ft and B form the basis for any fermion operator on the lattice, we 
summarize the basic properties of Q and B. 

1 The continuum limit in this paper stands for the so-called "naive" continuum limit with a — > 0, and 
the lattice size is gradually extended to infinity for any finite a in the process of taking the limit a — * 0. 
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2.1 Operators Q and B and Brillouin zone 

For a square lattice, for which we work in this paper, one can explicitly show that the 
simplest lattice fermion action 

S = J^fji^ (2.6) 
is invariant under the transformation [8] 

i/>' = Tip, $ = ^T- 1 (2.7) 
where T stands for any one of the following 16 operators 

1, TxT 2 , TxT 3 , T\T 4 , T 2 T 3 , T 2 T 4 , T 3 T 4 , TiT 2 T 3 T 4 , (2.8) 

and 

T 1; T 2 , T 3 , T 4 , T ± T 2 T 3 , T 2 T 3 T 4 , T 3 T 4 T 1; T 4 TxT 2 . (2.9) 
The operators T M are defined by 

T M = 7 M 7 5 exp(m^/a) (2.10) 

and satisfy the relation 

T li T„ + T,,T li = 26 lu , (2.11) 

with Tj = T M = T^ 1 for anti-hermitian 7^. We denote the 16 operators by 7^, n = ~ 15, 
in the following with % — 1. By recalling that the operator T M adds the momentum 7r/a 
to the fermion momentum k^, we cover the entire Brillouin zone 



— < kn < — (2.12) 
2a 2a 



by the operation (2.7) starting with the free fermion defined in 

-— <K<— . (2.13) 
2a 2a v ; 

The operators in (2.8) commute with 75, whereas those in (2.9) anti-commute with 75 
and thus change the sign of chiral charge, reproducing the 15 species doublers for (2.6) 
with correct chiral charge assignment; Sn=o( — 1)™75 = 0- 

One may define the near continuum configurations by the momentum k^ carried by 
the fermion 

- T^e < A;„ < -^e (2.14) 
2a ~ M ~ 2a v ' 

or 

-^e<ak^<^e (2.15) 

for sufficiently small a and e combined with the operation T n in (2.8) and (2.9). To identify 
each species doubler clearly in the near continuum configurations, we also keep r/a and 
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m /a finite for a — > small [8], and the gauge fields are assumed to be sufficiently smooth. 
For these configurations, we can approximate the operator D w by 

D w = ifl + M n + 0(e 2 ) + 0(a(<?^) 2 ) (2.16) 

for each species doubler, where the mass parameters M n stand for Mq = — ^ and one of 



2r m 4r m 

, 4,-1); , 6,1 

a a a a 

5:-=!, (4,-1); (1,1) (2,17) 

a a a a 

for n — 1 ~ 15. Here we denoted (multiplicity, chiral charge) in the bracket for species 
doublers. In (2.16) we used the relation valid for the configurations (2.15), for example, 

D w e ik * = J2 D w(x, y y k y 

y 

= [ j2Y^^ + -Y,( 1 - cosa K)-—V kx 

. . CL CL . . CL 

= [ 7 ^(l + 0(e 2 )) + -0(e 2 )-^]e^ (2.18) 

(X (X 

for vanishing gauge fields. 

For the near continuum configurations, we thus have from (2.3) 

D<$ +1) = i{]/)) 2k+l + M T f +1 ) + 0(e 2 ) (2.19) 

where the mass parameters M^ k+V> stand for 

M (2*+l) = _ ( ^0 )2fc+ l (2>2Q) 



a 



and one of 



(-) 2fc+1 - (^) 2fc+1 , (4,-1); (-) 2fc+1 - (^) 2fc +!, (6,1) 
(x (x a a 

(-) 2fe+1 - (^) 2fc+1 , (4,-1); (-) 2fe+1 - A 2 *^ 1 , (1,1) (2.21) 
a a a a 

for n = 1 ~ 15, in the same notation as in (2.17). 

To avoid the appearance of species doublers in 75-D, we choose M^ k+V> < and all 
other mass parameters M^ k+1 ^ > 0, n^O, namely 

< m < 2r. (2.22) 

The choice 

2m 2fc+1 = 1 (2.23) 
normalizes properly the Dirac operator H^k+i) in (2.2) 

H i2 k+i) ^ Hl^ Pf h+1 + iSl^ fl) 2{2k+1) (2.24) 
in the near continuum configurations for all \M n \ — > large [1]. 



3 Evaluation of the lattice Jacobian 

For an operator 0(x,y) defined on the lattice, one may define 

O mn = <f>m( X ) ( X i Vj&niv), (3.1) 



x,y 



and the trace 



TrO = J2°nn 

n 

n x,y 

= Y,CE<t>*n(x)0(x,y)<f>n(y))- (3-2) 

x n,y 

The local version of the trace (or anomaly) is then defined by 

trO(x, x) = ]T <P* n (x)0(x, y)4> n {y). (3.3) 

n,y 

For the operator of our interest, we have 
trT 5 (x) 

= tr[ l5 - ( l5 aD) 2k+1 ] 
= -tr( l5 aD) 2k+1 

^-^ 75[1+<+1 V(< +l Wr i) ] 

= -tr\[D% k+1) , 1 ] (3.4) 

-_!ffr ^ r -ikx r -l D (2k+1) 1 r ikx 

where we used the plane wave basis defined in the domain (2.13) combined with the 
operation T n . In this calculation, we repeatedly used the relation 

tr 75 = (3.5) 

which is expected to be valid in lattice theory. We also used a short hand notation 

Oe ikx = Y,0(x,y)e lky . (3.6) 

y 

There are various ways to evaluate the above trace (3.4). We evaluate the trace (3.4) by 
following the procedure used for the overlap Dirac operator in Refs. [9] [10]. Some of the 
basic papers of the lattice anomaly calculation are found in [11]- [14]. In this section we 
simplify the expression of the Jacobian, and its explicit evaluation is presented in the next 
section. 
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3.1 General analysis of the trace 

Our starting formula is (by using the momentum domain (2.12)) 

-m^^^^F^W^ (3 ' 7) 



with 



bf^(p) = (a^)D^%) 



= i[i^2>fwi Pli + a(?] 2fe+1 + [r£(l - cosp M ) + aB] 2k+1 

-(m ) 2fc+1 (3.8) 
and we defined the integration variable 

Vp, = ak^. (3.9) 

We used the definitions 

e~ ikx a $e ikx h(x) = [i ^ 7 M sin ak^ + cu$]h(x) (3.10) 

and 

e- ikx aBe ikx h(x) = [r ]T(1 - cos ak^) + aB]h(x). (3.11) 

for a sufficiently smooth function h(x). In the following we often omit writing h(x). 
Consequently, 

D {2k+1) (k,) = e~ ikx D {2k+1) e ikx 

= i[iY, + $] 2k+1 + [- £(1 - cos ak,) + B] 2k+1 

n a a n 

_ ( ^o )2fc+ i (312) 
a 

and Dy^ +1 \p) is defined by setting k^ = p^/a in D ( ^ r k+1 \k fl ). 

In the continuum limit a — > with p M = a/c M kept fixed, the operator approaches 

q = J2Y (cos ak^dp + ig cos ak^A^) + 0(a) = ]T 7" cos p MJ D M + 0(a) (3.13) 

and the leading term of B is known to be [11] 

B = -ir^smak^D^ + 0(a) = -irJ2 sin Pv D n + °(a) (3.14) 

with the covariant derivative defined by 

D^dp + igAp. (3.15) 
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Note that the "conventional naive continuum limit" is defined by a — > with kept 
fixed, instead of PfJi = ak^ being kept fixed as in the above limit. 
In the denominator of (3.7), one has a factor 

= [(^ 7 M sinp M + a<?) 2 ] 2fe+1 
+ {[ r J2(l-cos Pfl ) + aB] 2k+1 - (m ) 2fe+1 } 2 
-#£ 7 M sinp M + a$] 2k+ \ [r^l - cosp M ) + aB] 2k+1 ] 

= {£(sin^ - aiC»f + ^ [Y,Y}[C»,C„]} 2k+1 
+{[r £(1 - cos^) + a£] 2fc+1 - (m ) 2fe+1 } 2 

7" sinp„ + a$] 2k+1 , [rJ2(l - cosp M ) + aB] 2k+1 }. (3.16) 

Note that the first two terms in this last expression commute with 75, while the last 
term anti-commutes with 75. The last term of (3.16) is the interference term: From the 
structure of the commutator, one can confirm that it consists of terms with a factor 

-ia 2 Y[C„,B] = -mV^ 7 ^cos^sin^F^ + 0(a 3 ) (3.17) 

and the 2k factors of [«I] M 7 M sinp M + o0\ and the 2k factors of [rXI M (l — cos Pfl ) + aB}. 
We also note that 

[C^C U ] =igce>8p ll <x>8p v F tlv + 0(a). (3.18) 
To simplify various expressions in the following, we define the variables 

c M = cos akn = cos PfJi , s M = sina/c M = sin PfJi (3.19) 

and 

£7*011^, s 2 = 5>„) 2 . (3.20) 
3.2 Contribution of mass terms 

We now examine the integrand of (3.7) with only the "mass terms" in the numerator 
retained 

V 75 {[r E(l - S) + aB] 2k ^ - (m ) 2 ^} 1 . (3.21) 

V(^ +) (P)) f ^ +) (P) 

The numarator contains no 7^'s. We expand the denominator in powers of the interfer- 
ence term, which contains an odd number of 7 M 's. By remembering that the rest of the 



denominator factor contains an even number of 7 M 's and thus commute with 75, the odd 
powers in the interference term in this expansion anti-commute with 75 and thus vanish 
after taking the trace. Only the even powers in the interference term could survive the 
trace operation with 75. Since the interference term is of order 0(a 2 ), only the zeroth 
order term and the second order term in the interference term are important in the limit 
a -> 0. 

We can see that the second order term in the interference vanishes. Since the second 
order term is already of order 0(a 4 ) because of (3.17), one can set a = in all the 
remaining 2k powers of [i §/+ a(jH\ and the 2k powers of [r 5^(1 — c M ) + aB]. Namely these 
terms are replaced by i ^ and r^A\ — c M ), respectively. Since these factors commute 
with each other, the interference term consists of a sum of terms of the structure 

(i 9)'ioVE^(i ^>E(! " c m)P (3-22) 

with I + m — 2k, if one uses (3.17): Thess terms are linear in 7^ if one uses the relations 

(i = , 2 , 

(i$fr"M = -2($K-sV- ( 3 -23) 

We thus have only two 7 M 's with order 0(a 4 ) in the numerator of (3.21), which vanishes 
after the trace with 75. Note that the terms which contain 7 M 's in the denominator is of 
order 0(a 2 ) as in (3.16), and thus these terms cannot be used to supply extra 7 M 's. 

We have thus established that only the zeroth order term in the interference survives 
in (3.21), namely, one can set the interference term to in the denominator. In this case, 
from the expression in (3.16) we see that the 7^ factors appear only in the combination 

jWn u }[C»A). (3.24) 

The second power in this factor is just sufficient to survive the trace operation and cancel 
1/a 4 in front of the integral. This means that we can set a = everywhere except in the 
prefactor in (3.24). The expression (3.21) is then replaced by 

WS 1 - c ^ k+1 - (™o) 2fc+1 }^= (3.25) 



where 



a 2 



+{[r 2(1 " c,)f k+l - (m ) 2fc+1 } 2 . (3.26) 



3.3 Contribution of kinetic term 

Similarly, we analyze the integrand of (3.7) with the "kinetic term" in the numerator 

V75W affi** 1 } \ ■ (3.27) 
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Since the numerator is now odd in powers of 7 M , only the odd powers of the interference 
term could survive. The third power of the interference is 0(a 6 ), and only the first power 
in the interference need to be analyzed. 
We first rewrite the numerator factor as 

[i §/+ a^] 2k+l (3.28) 

= {5>„ - aiC») 2 + ^kf,-f][C„C v ]} k [i ^+ afl 

which shows that we have only one 7 M which is not multiplied by a. As we have already 
explained, the interference term in the denominator 

- ^+ aqf k+ \ [r£(l - c,) + aB] 2k+1 ] (3.29) 

is written as a sum of terms with a single commutator 

- ia 2 Y[C^ B] = -ia 2 gr ]T 7%*^ + 0(a 3 ) (3.30) 

multiplied by the 2k factors of [i ^+ aQ\ and the 2k factors of [r — c ^) + a B] ■ In such 
a term, if one exchanges the order of [i $f + aQ\ and [rX) M (l _ c m) + a -^]> one generates 
another commutator as in (3.30). We then have a factor a 4 and thus we can set a = 
in all other terms in the integrand. We recognize that such a term contains 7 M only in 
the combination with q/ and the two factors of the above commutator (3.30). From this 
combination together with / in the numerator (3.28), we cannot form a non- vanishing 
contraction with the antisymmetric e' iual3 tensor. 

This means that we can write the interference term as a sum of terms of the structure 

[i ^+ a$] l (-i)a 2 gr £ -fc^F^i ^+ aQ\ m 

x[rJ2(l-c,) + aB] 2k (3.31) 

where I + m — 2k. If both of I and m are even, we can use 

[i a$} 2 = - aiC,) 2 + £[f, C v \ (3.32) 

and 7 M always appears in the combination except for the numerator term (3.28) which 
contains and the commutator term (3.30), which contains a 2 7 M . Such a combination 
could give rise to a non- vanishing result. On the other hand, if both of I and m are odd, 
one has to deal with a left-over term 

[i t/+ aQ\{-i)a 2 gr ^ [i §/+ a$\ 

= [i t/+ a<y](-2)(-i)a 2 gr^c^s v F^[is ^ + aC^] 

fl,U 
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+[i f/+ a@}(-i)a 3 gr £ ^^c^c^D^F^ 



jXV 

H,u,a 

-{5>„ - aiC,) 2 + i 7 ^ 7 1[^,a]}(-i)aVE7%^- 



(3.33) 



where we used 7^7° + 7^7^ = —2r]^ a . The first term of this equation contains the factor 
a 2 [i y+ a(?], which should be replaced by a 2 i and should be combined with the factor % §/ 
in the numerator factor (together with a 2 [ 7 M , 7"][Cyt, C u ] from other factors) to obtain a 
possible non-zero result. But such a term cannot make a no n- vanishing contraction with 
e^ ualS . The second term is of orfer 0(a 3 ) and contains 3 7 M 's. If this term is combined with 
[i q/-\-a(/!!\ in the numerator, it becomes of order 0(a 5 ) due to (3.32). If it is combined with 
the drivative operator such as aC^ in (3.16), when commuting with other denominator 
factors, it becomes 0(a A ); such a term contains ^ = —s 2 if combined with / in the 
numerator and vanishes after trace with 75. Thus we can set the first and second terms 
to in the above equation (3.33). 

Only the last term in (3.33) can survive: Among those surviving terms, the even / 
terms and odd / terms cancel pairwise except one term in the interference term. 

By this way, we can write the total interference term as 



a 2 

-(2Hl)iflVE7"c M ^{« 2 + ^T[/j1v^]} fc 
x[rE(!-^)] 2fe (3.34) 

where the factor 2k + 1 comes from the 2k + 1 powers of [r — c M ) + aB\. We also set 
a = in all the terms without 7^ since this does not influence the surviving terms. We 
can also set a = in the numerator factor except for the combination a 2 [7 M , 7"]. Note that 
the order of 7^ and [7^, 7"] can be changed freely in the expansion of the denominator in 
powers of a 2 , since the surviving terms are contracted with 75 to give rise to e^ vafi ' . 

To summarize this tedious analysis, we can write the integrand with the "kinetic" 
term (3.27) as 

±tr l5 {(2k + l)rfa£)[7* 7%c^(£ i) 

x{ S 2 + ^^[ 7 ^7lc,c^]} 2fc [r5:(l- CM )] 2fc }-^ (3.35) 

where F^) is defined in (3.26). In writing this final expression we used the following 
sequence of rewriting 

' /'y ' I^CfiSvFfiu ^ ' ^ ] 7 / y^c^s u s a F^ l/ 
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= "^E^tIv^-E^O (3.36) 

Namely, only the term with two 7^'s contributes to the final result, and the odd term in 
s v after integration over the momentum vanishes. In the last step, we used the lattice 
hypercubic symmetry by taking into account the contraction with the e^ vafi symbol later. 



4 Formula for the chiral anomaly and parameter in- 
dependence 

The basic formula for the chiral anomaly is given by (3.7), (3.25) and (3.35). The next 
step is to expand the integrand in the powers of 

[ig^yCi^ -i JiCfiCvF^ and retain only the terms which contain the second power of this 
factor. We then combine the expansion with the formula 

tr^iiigh^^c^F^} 2 = (n c a )gHre^F» v F a p (4.1) 

4 a=l 

where e 1234 = 1. We thus write only the coefficients of the factor 

gHre^F, u F aP (4.2) 

in the following. 

The contribution from the "mass terms" (3.25) is given by 

1 d*p ,A , (2k + l)M (2fc+1) [3(2fc + l)(sr k - Ak{s^H] 

" 16 7- f (2n)^U a) H 5 / 2 [ } 

where 

H ee (sY k+1 + Mf 2k+1) 

M (2fc+1) ee [r ]T(1 - c M )] 2fc+1 - mf + \ (4.4) 
The contribution from the "kinetic term" (3.35) is written as 

x{«( S 2 ) 2 *- 1 ff-3(2* : +l)( s 2 )«} 7 ^ 7 ;. (4.5) 
Thus the total contribution is given by 
2A; + 1 /■¥ d A p 4 



'2fc+l — 



16 



/ ! 7|^(n^){M( 2fe+1 )-r(E^)[rE(l-^)] 2fe } 

J -2 a=l 8 C P H 



x{3(2k + l)(sr k -^(sr k ~ 1 H} 1 ^- 2 . (4.6) 
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4.1 Parameter independence 

To analyze the parameter independence of the coefficient of the chiral anomaly, we follow 
the procedure in Refs.[9] [10]. We first rewrite the integral in the domain — | < < ^ 
to the integral in the domain — | < p M < | by using the variables = sinp M as 

where the symbol takes care of the 16 (would-be) species doublers 

e„ = (±,±,±,±). (4.8) 

The following formula is also valid 

c, = ^{l-slfl\ (4.9) 

In this new notation, we have 

H{s) = (sr k+1 + Mf 2k+l) 

M (2k+1) = [r E(l " - S J) 1/2 )] 2fe+1 - mf^ (4.10) 

and we evaluate 

{2k + l)H(s) + \H{8)W E s^H{s)-^ 

= {2k + l)H(s) + iiJ( S ) 5 / 2 £ s^H(s)-^ 

= -{2k + l)M (2fc+1) {-[r E(l - " ^) 1/2 )] 2fc+1 + mf +1 

+[rEa - ^(i - ^r-E^a - sy- 1 / 2 } 

= (2fc + l)M (2fc+ i) (4.11) 
x{M (2fc+1) - [r£(l-£,(l " 4) 1/2 )] 2fe r£^(l - S J)-V2 } 

which is a generalization of the identity discussed in [10]. By using these relations one 
can prove 

1 d {s 2 Y k 

2^+T [ ^ s ^ + 4]( ^F 2) 



„2\4fc 



H{sf/ 2 

+{M (2fc+1) - [rE(l " e»(l - <) 1/2 )] 2 VE^(1 " sl)^ 2 } 



„2\4fc 



x5M (2fc+1) (-^^) (4.12) 
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and 



1 d (s 2 ) 2k - 1 



f 1 

2\2k-l 



His) 3 / 2 

+{M (2fc+1) - [rQl - 6,(1 - ^) 1/2 )] 2fe r£^(l - slV 2 } 



„2\2fc-l 



x3M (W^^)- ( 4 - 13 ) 



One can then show that 



dl 2k+l 2k + 1 J2(l[e ) fl 



dm 2k+1 16 ^ V /-I (2tt) 4 
x{3( 2A; + 1)1 W " 51^(2*+!) - r(E E(l " ^] 2 lM (2fe+1) ^-] 

( 2\2fc-l „2 _ / 2\2fc-l 

^M 1 ^- - Wtt+u - r (E ^ B 1 - c M)] 2fc ]M (2fe+1) ^ 7 ^]} 



2A; + 1 /TT x /-i 



zoi 



16 e ^ e »U(2*y 

Ak d .As 2 ™- 1 

+ 



s 



Similarly, we can show the relations by using (4.11) 

1 ^ _d_ , >E M (l-c M )] 2fc+ V) 4 \ 
2^ + 1^^ ^) 5/2 



{[r£„(l - c,)] 2fc+1 - [r - c,)] 2k rJ2,^sl(l - S 2 )-V 2 }( s 2 ) 



-ifc 



and 



H{sf' 2 

+{M (2fe+1) - [r £(1 - cgpr£ S ^(l - sj)" 1 /*} 



1 ^ _d_ , , u [rE,(l-c,)] 2fc+1 ( g 2 ) 2fc -\ 
2^ + 1^"^ +U #(s) 3 / 2 ] 

{VY^AA - c,)? k+1 ~ [rS„(l - ^)] 2fc rE,e, S 2 (l - s 2 )- 1 / 2 }( g 2 ) 2fc - 1 

H(sf/ 2 
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^UiA^pki, (4 , 6) 

where c M = e M (l — s^) 1//2 is understood. Using these relations, we have 

4, 



r 2fc+i dW 2fc+ 1 /-i rf 4 s 

9r 2fc +! 16 ^J- 1 } ^ J-i (2nY 



e M =± 

4fc 



v ro,o, . iu i!^ - c ^ )]2fc+1 - [rE ^ (l - - ^)- i/2 }( g2 ) 

-[M (2fe+1) - [rE(i - gi 2 VE^(i - 4r 1/2 ] 

x5M (2fc+1) ( #557/2 )} 

d , r {[r gg(l - c,)? k+1 ~ [r S„(l - c M )] 2fc r ^ 2 (1 - ^-^H* 2 ) 2 ^ 1 

1 ^( S )3/2 
2fc+l/„2\2fc-l 



[rE,(l-c,)] 2k+ \s 
H(sfl 2 



x3M (2fc+1) ( L ^ ^ s/9 v ' )}} 



2A; + 1 ^ ,„ s /-i d 4 s 



i6 J^n^u, 



7T 



The integrand becomes singular only if the following two relations simultaneously hold 
s 2 = 0, [rQl - 6,(1 - S J) 1/2 )] 2fc+1 - m 2fe+1 = (4.18) 

namely, only when m /r = 0, 2, 4, 6, 8. We are working in the physical region 

0<m <2r (4.19) 
and thus the above integrals are regular, and we have from (4.14) and (4.17) 

after partial integration. It can be confirmed that boundary terms at s M = ±1 give 
vanishing contributions after a summation over J2 e =±- This shows that the coefficient 
of the anomaly is stable under a smooth variation of the parameters r and itlq, which is 
expected for a topological quantity such as the chiral anomaly. 
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4.2 Explicit evaluation of the chiral anomaly 

Since the coefficient of the anomaly is independent of the parameters r and mo, we evaluate 
the anomaly in the limit where both of r and mo go to 0. To be precise we introduce an 
auxiliary parameter a, and take a limit a — > with both of 

T ~, (4.21) 
a a 

kept fixed in the physical region (4.19). The parameter a plays the role of an effective 
lattice spacing, though our formulas (4.3) and (4.5) are derived in the limit of the lattice 
spacing a = 0. 

There are various ways to evaluate the coefficient of the anomaly , and one of these 
methods is given in [10] in the analysis of the overlap operator with k = 0. We present 
a calculation which reveals a close connnection with the naive continuum limit. We first 
observe that the contribution of the "kinetic" term (4.5), which arises from the interference 
term in the denominator, vanishes in the above limit (4.21) To show this we examine 

x{4A;( S 2 ) 2fc - 1 iJ - 3(2k + l)( s 2 ) 4fc }-^. (4.22) 

where c M = e M (l — s 2 ) 1 / 2 is understood. We define the integration domain 

- e < s M < e (4.23) 

for all n with a sufficiently small but finite e. Since s 2 > and the denominator of the 
integrand is regular for the integration domain outside the above domain, the integral 
outside the domain (4.23) vanishes in the limit a — > 0. Note that the denominator of 

J2/3^ does not cause any divergence in the integral (4.22). In fact one can even take 
e — > in such a manner that 

a/e l -> (4.24) 

for a suitable fixed positive integer I. This is because the integral outside the domain 
(4.23) vanishes at least linearly in a, and thus one can let e — >• simultaneously with the 
above constraint, where the denominator e l stands for the possible infrared singularity in 
this calculational procedure. 

We thus examine the remaining integral 



^{Ak{s 2 f k ^H - 3(2k + l)( s 2 ) 4fc }-^. (4.25) 



This property has been used in the treatment of the overlap operator in Refs.[7] and [15]. 
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Since |c M | ~ 1 in the above integral, we can ignore the variation of c M . If one rescales the 
integration variable s M = as'^ and defines 

r' — r/a, m' Q = m Q /a (4.26) 

the above integral is written as 

~ £ 01 <*) /"; £Lv<L %>VTV - 

e M =± n J - e / a \ Z7T ) p C /3 /i 

x{4A;(( S ') 2 ) 2 ^ 1 ^ - 3(2A; + l)^)*)*}^ (4-27) 

where H(s') is parametrized by r' and m' , which are kept fixed in the limit a — > 0. 
Note that the factor (X)/? ^) i n the numerator, which is written in terms of the original 
variables, is 0(e 2 ). The above integral is convergent in this limit a — > and of order 
0(e 2 ), and thus it can be made arbitrarily small. We can even make it vanish precisely by 
taking the limit (4.24). We can thus ignore the contribution of the "kinetic" term, which 
arises from the interference term in the denominator, in the above limit (4.21) . 

We now come to the main contribution of the "mass terms" in (4.3). It turns out to 
be more convenient to go back to (3.25), which gives rise to (4.3). If one uses the notation 
of (3.4) instead of (3.7), (3.25) is written as 

E(^ L )(- 1 )4 r, S 75 {[rE(l±^)f W -N 2W } , ) -- (4.28) 
n=o 2 a 4 7-./2 (2tt) 4 m ^/F {k) (n,p^ 

where 

F {k) (n,p„) = {(t^+Y,aYc,D,) 2 } 2k+1 

+{[r]T(l± C g] 2fe + 1 - (m ) 2fc+1 } 2 . (4.29) 
/* 

The summation runs over the 16 would-be species doublers and the factor r^ M (l ± c M ) 
arises from each momentum domain. For later convenience, we modified the denominator 
factor in (3.25) by replacing s 2 + ig^-[y^ : l^c^F^ with (i ^+ J2 fJ ,0' / y^c IM D fl ) 2 , but this 
does not change the result as was explained in detail in the passage from (3.21) to (3.25) 
in Section 3.2. 

In the present integral, we can also define the domain 

- e < < e (4.30) 

for arbitrarily small but finite e. One can again confirm that the integral outside this 
domain vanishes at least linearly in a for a — > 0. This is because we retain the second 
order term in a 2 [7 M ,7 iy ]F^ to survive the trace with 75, and this cancels the factor 1/a 4 
in front of the integral. The resulting integral is finite outside the above domain, and it 
vanishes at least linearly in a in the limit (4.21). We can also let e — > as in (4.24). 
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We thus examine (4.28) inside the domain (4.30) 

U 2 ,( l) a' L (2n)^ Jf^S 

= l^ )( - ir r « tni Mi^w!!: (4 . 31) 

n=o 2 ^/a(2vr) 4 ^F {k) (n,ak^) 

For sufficiently small e, we have 

[r£(l ± cosa£g] 2fe+1 /a 2fe+1 " (m ) 2fe+1 /a 2fc+1 = M< 2fe+1 ) + 0(e 2 ) 

F (fc) (n, aA; M )/a 2 ( 2fc+1 ) = (i ^+ + (M< 2fc+1 >) 2 + 0(e 2 ) (4.32) 

where the mass parameters Mff k+1 ^ are defined in (2.20) and (2.21). The above integral 
in the limit a — > with (4.24) approaches 

2 n=0 ( 27r ) ^([i ^+ ^)]2)2fe+l + (M^ fc+1) ) 2 

1 r°° d 4 k 1 

/ 7773775" 



2 7-oo (2tt) 4 v /((^ +j ^)2/M2)2fc+i + 1 

~E(-l) B *r/ 7^75 , (4.33) 

by recalling M < 0, and M 2 = [(M^ 2fe+1) ) 2 ] 1 ^ 2fc+1 ). The sum of integrals in (4.33) gives 
rise to the anomaly for all M 2 — > 00 in the final stage: [15] 

roo J7) 2 n 2 

lim tr / ^—je~ ikx -f 5 f(^--)e ikx = -*—tre! a ' af> *F U1/ F a0 . (4.34) 

m^oo i_oo (2vr) 4 ,5M M 2; 32tt 2 ^ 1 K J 

Here we defined 

fix) = - 1 (4.35) 



which satisfies 

/(0) = 1, /(oo)=0, 

/'(a:)x| x= o = /'(^xl^oo = 0. (4.36) 

The left-hand side of (4.34) is known to be independent of the choice of f(x) which satisfies 
the mild condition (4.36) [6]. 

By combining (1.15), (3.4) and (4.34), we recover the Atiyah-Singer index theorem ( 
in continuum i? 4 space) [16] [17] 

n + -n„ = j dx^tre^F^F al3 (4.37) 

in the (naive) continuum limit for any operator in (1.1), if one follows the construction of 
l5 D in [1]. 
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5 Discussion 



We have shown in an explicit and elementary manner that a new class of lattice Dirac 
operators proposed in [1] satisfy the correct anomaly relation. We have in particular 
shown that the anomaly coefficient is independent of a small variation in the parameters 
r and mo, which characterize these Dirac operators. This is in agreement with the more 
general but formal analysis in [1]. 

The Dirac operator in [1] is defined in a somewhat indirect manner as 

det#=(dettf( 2fe+1) ) 1/(2fc+1) (5.1) 

as in (A. 9) in Appendix. This definition is sufficient for the non-perturbative analysis, 
and as we have shown in this paper, it is also sufficient to evaluate the chiral anomaly 
explicitly. However, the perturbative treatment of these general class of Dirac operators 
is not well understood yet. It would therefore be interesting to extend the analyses in 
[14] [18], for example, to the present general class of operators. 

As for the locality issue of the present class of operators, the fact that the anomaly 
calculation in a naive continuum limit makes sense suggests that there is a certain range 
of gauge field configurations which make these operators local. The fact that we take 
a 2k+lth root of H^k+i) in (5.1) by itself may not spoil much of the locality, since the 
eigenfunctions in (A. 17) are defined in terms of H^k+i) and thus they may reflect the 
locality properties of H^k+i) , which may not differ qualitatively from those of the overlap 
operator [19] [20]. In any case, a direct analysis of this locality issue is left as an important 
problem. 

Our construction of the Dirac operator (5.1) on the basis of the defining algebraic re- 
lation (1.1) suggests that we obtain better chiral properties if one increases the parameter 
k. An intuitive argument for this expectation is that the right-hand side 2a 2k+1 ( , ~f 5 D) 2k+ ' 2 
of the algebra (1.1), which breaks chiral symmetry, becomes more irrelevant for larger k in 
the sense of renormalization group. At the same time, however, our construction requires 
a larger lattice for larger k since the basic operator appearing in our construction spreads 
over far apart lattice points for large k. To maintain the locality of the Dirac operator, 
we need to take a smaller lattice spacing for larger k. 

As for the chiral fermions[21] [22], the present calculation of anomaly is readily extended 
to the evaluation of the fermion number anomaly of chiral theory and also to the so-called 
covariant form of non-Abelian anomalies [6] in the continuum limit. But the construction 
of chiral fermion theory at finite lattice spacing is a challenging un-solved problem not 
only in our general Dirac operators but also in the original overlap operator [3]. 

T-W. Chiu has recently informed us that a numerical study of some of basic prop- 
erties of the operator with k — 1, such as index theorem, chiral anomaly and the 
propagator, is in progress [24]. 

One of us (KF) thanks Ting-Wai Chiu for numerous helpful discussions from the very 
beginning of the present investigation. 
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A Basic construction of general Dirac operators 

We start with (1.1) written in the form[l] 

H l5 + l5 H = 2H 2k+2 (A.l) 

or equivalently 

T 5 H + HT 5 = (A.2) 
where H = a^D and T 5 = 75 — H 2k+1 . This algebraic relation implies that 

75 # 2 = [75# + H l5 ]H - H[ lb H + H l5 ] + H 2 l5 = H 2 lb . (A.3) 

Namely, the algebraic relation (A.l) is equivalent to the two relations 

H 2k+1 lb + lb H 2k+1 = 2H 2{ - 2k+1 \ 

l5 H 2 - H 2 l5 = 0. (A.4) 
If one defines H^k+i) = H 2k+1 , the first relation of (A.4) becomes 

H(2k+i)l5 + 75#(2fc+i) = 2H 2 2k+l) (A. 5) 

with T 5 = 75 — H(2k+i), which is identical to the conventional Ginsparg- Wilson relation 
with k = in (1.1). We utilize this property to construct a solution to (A.l). 

The physical condition for the operator H in (A.l) in the near continuum configuration 

is 

H~ l5 aip + 75(75^ $>? k+2 (A.6) 

where the first term stands for the leading term in chiral symmetric terms, and the second 
term stands for the leading term in chiral symmetry breaking terms. Thus H^k+i) should 
satisfy 

H {2k+ i) * [ 75 « V> + 75(750* P) 2k+2 ] 2k+1 

~ (75^^) 2fc+1 + 7 5(750^) 2(2/£+1) (A.7) 

as can be confirmed by noting 75 p+ p^y 5 = 0. Here only the leading terms in chiral 
symmetric and chiral symmetry breaking terms, respectively, are written. 
One can thus construct a solution for H^k+i) by 



H, 



, a+1 , = ^[l + 7 5g fT" I — ] (A.8) 

V% H w 

in terms of the hermitian H^ +V> = 7 5 _D^ fc+1 " ) = (iJ^ fe+1 ^)t. The operator D^ k+1 ^ is 
defined in (2.3). The physical condition (A.7) is satisfied by (2.19), as was noted in the 
body of the text. 

We now discuss how to reconstruct H, which satisfies (A.l), from H^k+i) defined 
above. The basic idea is to define in the representation where H^k+i) is diagonal 

H = (H (2k+1) )^ 2k+ ^ (A.9) 
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in such a manner that H thus obtained satisfies the second constraint in (A. 4). For this 
purpose, we first recall the essence of the general representation of the algebra (A.l)[l]. 

If one defines the eigenvalue problem 

#(2fc+l)0n = (aA n ) 2fc+1 0n, {4>nAn) = 1 (A.10) 

one can classify the eigenstates into the 3 classes: 

(i) n± ("zero modes"), 

H(2k+l)<t>n = 0, 7 5 0n = ±<f>n, (A- 11 ) 

(ii) N± ("highest states"), 

H(2k+i)<Pn = ±0n, 750n = ±0n, respectively, (A. 12) 

(iii) "paired states" with < |(aA„) 2fe+1 | < 1, 

H (2k+l) ct> n = (aA n ) 2fc+1 0n, ^ ( 2 fe+ i)(r 5 0„) = -(aA n ) 2fe+1 (r 5 n ). (A.13) 

where 

r 5 = 75 -// {2fe+1) . (A. 14) 
Note that r 5 (r 5 n ) oc <f) n for < |(aA n ) 2fe+1 | < 1. 

We have a chirality sum rule [23] 

n + + N + =n_ + N^ (A. 15) 

where N± stand for the number of "highest states" in the classification (ii) above. 

All the states <p n with < |(aA„) 2fc+1 | < 1, which appear pairwise with (aA„) 2fe+1 = 
±|(aA„) 2fe+1 |, can be normalized to satisfy the relations 

T 5 n = [1 - (aA n ) 2(2fc+1) ] V V-n, 

750n = (aA^^Vn + fl-W 2 ^ 1 )] 1 ^^, (A.16) 
where 0_„ stands for the eigenstate with an eigenvalue opposite to that of <p n . 

We can define the solution H of (A.l) operationally by 

H<i) n = a\ n (f) n (A. 17) 

by using the same set of eigenfunctions and (the real 2k + 1th roots of) eigenvalues 

{0„}, {a\ n } (A.18) 

as for Hpk+i) in (A. 10). Note that the operator T 5 = 75 — H^k+i) — 75 — H 2k+1 , which 
reverses the signature of eigenvalues of "paired states" and defines the index, is identical 
to (A.2) and (A.5). 
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We can confirm the second constraint in (A. 4) and the defining algebraic relation (A. 2) 
for any n in (A. 17) by using (A. 16), 

[# 2 75 - 7 5# 2 ]0n = 

[r 5 # + #r 5 ]0 n = o. (A.i9) 

The general representation of the algebra (A.l) is obtained from the standard repre- 
sentation, which is defined by H and {<j) n } in (A. 17), and 75 in (A. 16), by applying a 
suitable unitary transformation. 
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